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Abstract. We investigate the initial- value problem for the relativistic Euler equations 
governing isothermal perfect fluid flows, and generalize an approach introduced by LeFloch 
and Shelukhin in the non-relativistic setting. We establish the existence of globally defined, 
bounded measurable, entropy solutions with arbitrary large amplitude. An earlier result 
by Smoller and Temple for the same system covered solutions with bounded variation that 
avoid the vacuum state. The new framework proposed here provides entropy solutions 
in a larger function space and allows for the mass density to vanish and the velocity 
field to approach the speed of light. The relativistic Euler equations become strongly 
degenerate in both regimes, as the conservative or the flux variables vanish or blow-up. 
Our proof is based on the method of compensated compactness for nonlinear systems of 
conservation laws (Tartar, DiPerna) and takes advantage of a scaling invariance property 
of the isothermal fluid equations. We also rely on properties of the fundamental kernel 
that generates the mathematical entropy and entropy flux pairs. This kernel exhibits 
certain singularities on the boundary of its support and we are led to analyze certain 
nonconservative products (after Dal Maso, LeFloch, and Murat) consisting of functions of 
bounded variation by measures. 
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1. Introduction and main result 

1.1. Purpose of this paper. The relativistic Euler equations describe the dynam- 
ics of a compressible fluid in the context of special relativity, i.e. a fluid evolving on 
the flat Minkowski spacetime. This system can be regarded as an approximation to 
the Euler-Einstein equations, which is valid in a small region of the spacetime and 
away from large matter concentrations. Under the assumption of plane-symmetry, 
the fluid unknowns consist of the co-moving mass density p > and the velocity 
field v € (— c, c), where c denotes the speed of light. Our main purpose in the 
present paper is the construction of weak solutions containing shock waves and 
having arbitrary large amplitude. 

The mathematical properties of this model were investigated in works by Taub 
(1957), Thompson (1986), Lichnerowicz (1993), and others. These equations are 
also important in computational physics and we refer to Marti and Miiller (2003) 
for an extensive review of the Riemann problem and the numerical methods in 
hydrodynamics in the context of special relativity. The mathematical analysis of 
the relativistic Euler equations for the isothermal fluids, considered in the present 
paper, has received much less attention in the literature; our emphasis is on fluids 
governed by the linear pressure law 

p(p) = k 2 p, 

where k > is the (local, constant) sound speed and by the principle of special 
relativity must be less than the speed of light denoted here by c. 

The Euler equations form a nonlinear system of partial differential equations 
of hyperbolic type. It is well-known that solutions even they are smooth initially 
will eventually become discontinuous and must be understood in the weak sense 
of distributions. Furthermore, for the sake of solutions it is necessary to constrain 
these weak solutions by certain entropy inequalities. 

The initial-valued problem for relativistic isothermal fluids was first studied by 
Smoller and Temple (1993). They established the existence of entropy solutions, 
under the assumption that the initial mass density is bounded, is bounded away 
from zero, and has bounded variation. Their result is based on the Glimm scheme 
and extends an earlier approach for the non- relativistic version (Nishida, 1968). 

In the present paper, we propose an alternative approach based on the method of 
compensated compactness for nonlinear conservation laws (Tartar, 1983) and, for 
the relativistic Euler equations, we provide a mathematical framework encompass- 
ing a large class of weak solutions. These solutions may take vacuum values p = 
and may have high- velocity approaching the speed of light. Solutions of this nature 
arise naturally in applications; for instance, a star is described by a compactly sup- 
ported mass density function. As can be checked easily, at points where p vanishes 
the Euler equations are highly degenerate and the conservative variables vanish 
identically while the velocity field is ill-defined. Another related singularity of the 
equations is obtained in the limit when the fluid velocity approaches the light speed 
and the wave speeds approach each other; hence, the model fails to be uniformly 
strictly hyperbolic. These features of the relativistic Euler equations for isother- 
mal fluids lead to particularly challenging mathematical questions, concerning the 
existence and the behavior of entropy solutions. 
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1.2. Relativistic fluid equations. We consider the following system of two con- 
servation laws 



(p + pc 2 )v\ iQ ( (p + pv 2 )c 2 ' = „ 



which, by setting e = 1/c and using the condition p(p) — k 2 p, reads 

„ fl + e i k 2 v 2 \ n fl + e 2 k 2 
d t [ 

(1.1) 



1 — e 2 v 2 J V 1 — e 2 v 



r fl + E 2 k 2 \ f k 2 +V 2 \ „ 

^{r^ pv ) +dx {r^ 2 ^ p ) =0 - 

Taking the formal limit e — 0, we obtain the non-relativistic version of these equa- 
tions: 

d t p + d x (pv) = 0, 
(1 ' 2) d t (pv) + d x ((k 2 + v 2 )p) =0. 

Our starting point is the recent work by Huang and Wang (2003) and LeFloch 
and Shelukhin (2005). The existence of entropy solutions to Q1.2p is known when the 
mass density is bounded and the velocity is unbounded. It was observed that the 
natural function space associated with the non-relativistic Euler equations allows 
for the velocity field to be unbounded. Our aim will be here to generalize to 
relativistic fluids the method and results by LeFloch and Shelukhin (2005), based 
on such "natural" estimates. 

The equations (II. 1|) form a nonlinear hyperbolic system of partial differential 
equations of the form 

(1.3) d t G + d x H = Q, d t H + d x F = 0, 

where 



1 + eW l + e 2 k 2 

G{P, V) = — Y1~P' H \P> V > = "I 2^2 P v > 

1 — e v 1 — e v 

k 2 + v 2 



F (P> V ) = T1.P- 

1 — £ V 

Since shock waves are known to arise even from smooth initial data (Pan and 
Smollcr, 2006), we need a concept of solutions that include discontinuous functions. 
We introduce the following notion of entropy solution. 

First of all, we will say that a Lipschitz continuous map (IX, 5F) is an entropy pair 
if every smooth solution of (|1.1[) satisfies the additional conservation law 

d t U(p,v)+d x ?(p,v) =0. 

However, weak solutions are required to satisfy the above conservation laws for 
convex functions U, but as inequalities only. Furthermore, the class of entropy 
functions is further restricted as we only consider weak entropies that is functions 
XL vanishing on the vacuum line p = 0. For instance, both pairs (G, H) and (iJ, F) 
in (|1.3p are (trivial) weak entropy pairs. We also set 

2e ._ . 
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Definition 1.1 (Notion of entropy solutions). A tame region is a set of the form 



for some constant M > 0. A pair of measurable and bounded functions Po,Vq : K — > 
R is a tame initial data if its range is included in a tame region. 

Given a tame initial data po,vo, a pair of measurable and bounded functions 
p, v : R+ x R — > R is called a tame entropy solution to the isothermal relativistic 
Euler equations (jTTTJ) associated with the initial data po,vo if the range of p,v is 
included in a tame region and, moreover, 



for every convex, smooth, weak entropy pair (XI, 5") of the isothermal relativistic 
Euler equations and every non-negative test-function 9 supported in [0, oo) x R. 

Observe that the inequality in the definition of a tame region allows the velocity 
to approach the light speed (normalized here to be when the mass density 
approaches zero. The concept of a tame region is quite natural, as it is equivalent 
to uniform bounds on the Riemann invariants (defined later in Section [2|). Clearly, 
the entropy pairs under consideration need not be globally Lipschitz continuous, 
but only so within any tame region, as those are the only regions of interest. 

One of our main results in this paper is the following existence theorem. 

Theorem 1.2 (Existence theory). Given any parameter value e £ (0, 1) and a tame 
initial data po,vo, the initial-value problem for the relativistic Euler equations for 
isothermal fluids admits a tame entropy solution p, v : R+ xl^R associated 

with Po,Vq. 

1.3. Main ideas for the proof. The approach proposed in the present paper 
will rely on the following observation, which was already pointed out in the non- 
relativistic setting. 

Lemma 1.3 (Linearity property). If(p,v) is a (weak, entropy) solution of the rel- 
ativistic Euler equations for isothermal fluid (jl.ip . then for every positive constant 
X, the function (Xp,v) is also a (weak, entropy) solution of the same equations. 

Our general strategy of proof follows, on one hand, DiPerna (1983), who obtained 
bounded solutions p > and »el for non-relativistic polytropic fluids satisfying 
p(p) = k 2 p 7 , with 7 > 1 and, on the other hand, LeFloch and Shelukhin (2005), 
who extended DiPerna's analysis to include isothermal fluids satisfying p(p) = k 2 p 
and observed that the velocity field v need not be bounded. 

The main difficulty for our analysis in this paper lies in the lack of uniform strict 
hyperbolicity of the Euler equations when the fluid velocity approaches the light 
speed. To deal with this problem we will proceed along the following lines: 

• Mathematical entropy pairs. Our first task will be to construct entropy 
pairs which amounts to solve a linear hyperbolic equation in the variable 
p, v. We will introduce the Riemann function 51 and the entropy kernel \ 
associated with this equation, so that the entropy pairs of interest can be 
expressed by an explicit formula in term of the kernel \. Contrary to the 
case of polytropic fluids (DiPerna, 1983) and in agreement with the case 
of (non-relativistic) isothermal fluids, an initial data for the entropy kernel 
must be imposed away from the vacuum, say on the line p = 1. 



T e (M) :={p,v/0<p< M, 1 - e|v| > [p/Mf'} 
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• Singularities of the entropy kernel. We will show that the function x is 
discontinuous along the boundary d% of its support, so that its first-order 
derivatives exhibit Dirac masses. In the case e = the kernel was given by 
an explicit formula. In contrast, when e > we need to derive uniform esti- 
mates on x and determine explicitly the traces of its first-order derivatives 
along d%. 

• A priori bounds. We will next derive a priori bounds on approximate 
solutions p h , v h generated by the Lax-Friedrichs scheme. Our bounds show 
that the mass density remains uniformly bounded and the velocity field 
satisfies the tame condition. 

• Reduction of the Young measure. In the final part of the proof we identify 
the structure of a Young measure v = v t , x associated with the sequence 
p ,v 1 . We analyze certain nonconservative products (Dal Maso, LeFloch, 
and Murat, 1995) consisting of functions of bounded variation by measures. 
The term of interest contains a key coefficient (denoted below by S(p)) 
which docs not vanish, provided we take advantage of the scaling invariance 
property in Lemma 1 1.31 



2. Basic properties of the model 

2.1. Wave speeds and Riemann invariants. Scaling properties of the equations 
(|l.ip will play an important role. Observe that the transformation 

v' = v/k, t' = kt, e' = ke 

allows one to reduce the system (jl.ip to the same system with k = 1. In view of 
the physical constraint < k < c between the sound speed and the light speed this 
amounts to impose c > 1. The limiting case c — > 1 corresponds to the special case 
where the sound speed and the light speed coincide. From now on, we suppose that 
k = 1 so that the Euler equations read 

l + e i v 2 \ „ / 1 + e 2 



Q t I i P)+9x[ 1 _ £ 2 V 2 P V ) = > 



The velocity is restricted to lie in the interval (— 1/e, 1/e); note that the conservative 
and flux variables blow-up when v — > ±l/e. The range of physical interest for e is 

< e < 1, 

the limiting case e = and e — 1 corresponding to the non-relativistic model (speed 
of light is infinite) and the scalar field model (the sound speed and the light speed 
coincide), respectively. 

Indeed, the system (|2.1[) in the limit e — > 1 converges to 

(2.2) v 7 v 7 

dt {ih^ pv ) +dx (\^ p ) =0 ' 
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which is simply equivalent to the linear wave equation. This is clear by introducing 
the unknowns a :— p and b :— j— 3 pv, so that 

dta + d x b = 0, 
9*6 + c^a = 0. 

The conservation laws (12. ip form a nonlinear hyperbolic system whose Jacobian 
matrix admits the two eigenvalues 

v — 1 v + 1 

Ai :— 7. — , A2 



1 — e 2 v ' 1 + e 2 v 

Clearly, the characteristic speeds are smooth functions in the closed interval v £ 
[—1/e, 1/e]. The corresponding eigenvectors are 

-1 1 V 



>'2 



1 — e 2 v 2 ' 1 + e 2 p 
1 1 1 



)■ 



. 1 — e 2 v 2 ' 1 + e 2 p/ 

The Riemann invariants u>, z, by definition, satisfy Vu> ■ ri = 0, Vz ■ r2 = 0, and 
are uniquely defined up to the composition by a one-to-one map: 

„ „ w — z 
w := u + R, R — , 

(2.3) ± 

z := u — R, u — , 

2 

where i?, u are functions of p, v. 



1 

hi 

1 , /l + ev\ le 2eu 



R = R{p);=-—\np, p=e^ + ^ R , 
1 + e z 



m = u (v) := — m 

V ; 2e \l-evJ 



The Riemann invariants provide a change of variables (fi,v) 1— > (u>, z), which will 
be often used. 

Clearly, the mapping v 1— > tt(w) is one-to-one from the bounded interval (—1/e, 1/e) 
onto the real line R. The mapping p 1— > R(p) is one-to-one from (0, 00) onto M. It 
is not difficult to check that, in terms of the variables w,z, the system (12. ip takes 
the diagonal form 

+ A2 d x w = 0, <9 t z + Ai d x z = 0. 

Observe that 

ll + e-(l-e)e £ ( lu+z ) 
M^)-— £l + E+{1 _ e)ee{w+z y 

(2-4) 

l l-e-(l + e)e^> 

»,Z = - 7 : — r = — Ai (— W, —Z). 

v ' ; e I - e + {I + e)e^ w + z ) n ' ' 
Sometimes, we will also use of the "modified" Riemann invariants defined as 

(2 _x ' Vl-ew/ 

^:= e -^pV(^)(i±^)- 1/(2E) . 
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Note that p = (WZ)^ 2 ^ 2 . 

Expressing now the physical variables p, v as functions of the Riemann invariants, 



p = R ((>-z)/2) = exp (l + £) 



and 



v = v(w + z) 



^ e e(w+z) _ ^ 



e e(w+z) _|_ l 



1 



e e{w+z) _|_ l 

we obtain the w- and z-derivatives of (p, v): 

1 1 



Pw = -pz 



tanh ( e 



w + z 



2 „,2 



= ^ 1 - e « 
2u„ 2 V 

1 1 + e 2 



Moreover, the derivatives of the Riemann invariants considered as functions of (p, v) 
are 

1 1 

Up = -z 



l + e 2 p 1 
1 



»'<; 



1 — e 2 v 2 

Finally, for latter use we express the Euler equations in the nonconservative 
variables (p,v). By setting 



D(G,H) 
D(p,v) 



/ 1 + e 4 v 2 2e 2 (l + e 2 )pv \ 

1 — e 2 v 2 (1 — e 2 v 2 ) 2 

(l + e 2 )v (l + e 2 )p(l + e 2 v 2 ) 



V i 



(l-e 2 v 2 ) 2 ) 



and 



D(H,F) 
D(p,v) 

we can rewrite (12.11) as 



I {i + e 2 )y (1 + £ >(1 + £ V) ^ 

l-e 2 v 2 (l-s 2 v 2 ) 2 

l + v 2 2(l + e 2 )pv 

V 1 - e 2 v 2 {l-e 2 v 2 ) 2 ) 



d t u + d x G(u) = 0, u 



with 



DG _ { D{G,H) \ - 1 D(H, F) 

~m ~ \ D(p,v) J D(p,v) 

( (l-e 2 )v (l + s 2 )p \ 

1 - e 4 v 2 1 - e A v 2 

(l-e 2 v 2 ) 2 (l-e 2 )v 



V (l + e 2 )p(l~e 4 v 2 ) l-e 4 v 2 J 
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2.2. Strict hyperbolicity fails for high-velocity fields. From the relation 



1 - e z v z 
A 2 - Ai = 2 > 



we deduce that: 



Lemma 2.1 (Hyperbolicity properties). The Euler equations for isothermal rela- 
tivistic fluids are strictly hyperbolic in the region \v\ < 1/e for all p > (i.e. even 
in the presence of vacuum singularities in the mass density field), but strict hyper- 
bolicity fails as v — > ±l/e (i.e. in the presence of light speed singularities in the 
velocity field). 

In contrast, for polytropic perfect fluids the Euler equations also fail to be strictly 
hyperbolic at the vacuum p = 0. 

The ^-derivatives of the eigenvalues considered as functions of (p, v) are 

1 - e 2 _ 1-e 2 

while Aip = X-2 P = 0. Their derivatives along the characteristic fields are 

1-e 2 1 

VAl ■ n = — 2 ~T\ T~Y2 > ' 

1 + e z p(l — e z v) z 

1-e 2 1 
VA 2 • r 2 = 9 2 2 > 0, 

1 + e z p(l + e z v) z 

while in terms of the Riemann invariants we have 

. _. . _ (l-£ 2 )(l~e 2 ^ 2 ) _ 

Aim — Mv v w + *\pPw — ; — ~2 — Aiz, 

2{l-e 2 v) 

[ ' , , x (l-e 2 )(l-e 2 v 2 ) _^ 

A2z — *2vVz + ^2pPz — ; " — "2 — A2m- 

2 (1 + e z v) 

Hence, we conclude: 

Lemma 2.2 (Genuine nonlinearity property). The Euler equations for isothermal 
relativistic fluids admit two genuinely nonlinear characteristic fields in the domain 
p > 0, \v\ < 1/e. However, the genuine nonlinearity property fails in the limit 
v -> ±l/e. 

In contrast, for polytropic perfect fluids the Euler equations also fail to be gen- 
uinely nonlinear at the vacuum. 

We point out that the wave speeds and Riemann invariants are smooth functions 
even as e — > 0: 

\i = (v-l)(l + e 2 v + 0(e 4 v 2 )), 
A 2 = (v + l)(l-e 2 v + 0(e 4 v 2 )), 

and 

(2.8) u(v) =v (l + jv 2 + 0(e 4 v 4 ] 

Note that (|2.7p is uniform in the whole interval v 6 [—1/e, 1/e], while the remainder 
in (|2.8|) blows-up when v approaches the light speed. 
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3. Derivation of additional conservation laws 

3.1. Entropy equation. From the equations and for smooth solutions we can 
derive additional conservation laws, which will play a central role in the existence 
theory. By definition, a pair of mathematical entropy U = U(p, v) and entropy-flux 
5F = 3(p, v) provides a conservation law satisfied by all smooth solutions of Ijl.ljl . 
The entropy pairs are determined by the compatibility conditions 

V^Tj = XjVU-rj, j = l,2. 

Expressing It, ? as functions of w, z and relying on the properties of the Riemann 
invariants, these conditions are equivalent to 

(3.1) J w = A2 U w , J z = \ilL z , 

and imply an equation satisfied by the entropy U = 1l(w, z) only: 

(X 1 U Z ) W = (X 2 U W ) Z . 

That is, XL satisfies 

(3.2) U wz + - X2 \ U w - — ^- U z = 0, 

which we will refer to as the entropy equation. 

Using the formulas (|2.6|) . the coefficients in (|3.2[) are found to be 

Xlw = (l-£ 2 )(l + £ 2 v) 

A2-A1 4(1 -e 2 v) 

X 2z _ (l-£ 2 )(l-£ 2 ^) 
A2-A1 4(l + £ 2 v) ' 

and thus (|3.2p becomes 

(3.3) U wz + b(w + z) VL W + a(w + z) U z = 0, 
where the coefficients depend on w + z only. We have set 

(l-£ 2 )(l + £ 2 v) 



a(v) := -- 



4(1 -e 2 v) 



(1 - e 2 )(l ~ e 2 v) 
4(1 + s z v) 

in which v — v(w + z), and a :— aov. b :— b o TJ, therefore 

= 1-g 2 l-£+(l + g)e £g 

(3 4) 4 + 

- l-£ 2 l + e+(l-e)e* 

K0 = — i- e +(l+ e)e * = 

The equation (|3.3[) is a linear hyperbolic equation with smooth coefficients, and 
its solutions are generated by the Riemann function 3l(w' , z'; iu, z), defined for each 
fixed (w, z) by the Goursat problem associated with the adjoint operator: 

Xv-z' - (b(w' + z') %) w , - (a(w' + z') %) z , = 0, 

Ji w > (w , z; w, z) = a{w' + z) 3l(w' , z; w, z) on the line z = z, 

(3.5) 

H z >{w, z ; 10, 2) = &(w + z ) ft(-u;, z'; u>, z) on the line u/ = io, 
3?(u>, z; w, z) = 1. 
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The Riemann function allows us to solve the general characteristic value problem 

Uwz + b(w + z) IL W + a(w + z) U z = g, 
(3.6) 11(w, z) — tp(w) on the line z = z , 

11(11)' , z) — ip( z ) on the line w = w' , 

where <p, ip are prescribed characteristic data and g is a given source. Indeed, we 
have 

U(w, z) :=— <p(w) !R(w, z'; w,z) + — ip(z) %(w' ', z; w, z) 
w 1 

-&(w",z';w,z) ip w (w")dw" 
+ ( Ua(w" + z') z'; w,z)-)- % w ,(w", z'\ w, z)\ <p(w") dw" 

Jul' ^ I J 

+ f ~%(w',z";w,z)ijj z (z")dz" 

J z' * 

+ J (jj>(w' + z")0i(w', z" ;w, z) — —0i z '(w', z" ;w, z)^j ip(z") dz" 
g(w",z") R(w", z"; w, z) dw" dz" . 



3.2. Non-relativistic limit. We assume first that e = 0. Formally, when s — > 
in (|3.3p we obtain 

+\k- u °) = °- 

The Riemann function associated with this equation was constructed by LeFloch 
and Shelukhin (2005): 

(3.7) %°(w', z'; w, z) = e ((«-«')-(*-»'))/4 f(( w _ w >) (z - z ')) j 

where the function f° = f°(m) is related to the Bessel function of order and can 
be characterized as the solution to the ordinary differential equation 

mf mm + f m + f/lQ = 0, 

/» = 1, /°(0) = -l/16. 

One can check that, for every fixed (w,z), the function (w',z') i— > ^R°(w',z';w,z) 
defined in (|3 . T[) is the unique solution of the Goursat problem 

$L'z' ~ 7 ' ~^z') = o, 

%P(w',z;w,z) = eS w - w "> /4 on the line z' = z, 

%,°(w, z'; w, z) = e~( z ~ z ^ 4 on the line w' = w. 

The following description of the entropy kernel was also established. 

Theorem 3.1 (Entropy kernel of the non-relativistic Euler equations). Consider 
the isothermal non-relativistic Euler equations (If .2|) (with k — 1 after normaliza- 
tion). Then the junction 



'ft°(0, 0; w, z) = e^-^)/ 4 f°(wz), wz<0, 
0, wz > 0, 
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is a fundamental solution of the entropy equation 

X° wz + \ (xi ~ X° z ) = ~2S W=Q <8> 5 Z=0 . 
It is a function with bounded variation and 

X° w (w, z) = e' z ' 4 ( - (sgn*) S w=0 + e w '\-f{wz) + z f? n (wz)) l wz<0 ), 

(3.8) 4 

X° z (w, z) = e w ' 4 ( - (sgn w) 5 Z=0 + e - z l\--f{wz) + w f m {wz)) l wz<0 ), 

where t wz< Q denotes the characteristic function of the set {wz < 0}. 

It is is natural to define the entropy kernel by imposing data on the line p = 1. 
(Observe that w — z — correspond to (p, v) — (1,0).) Hence, the mathematical 
entropies of (jl.2p that vanish on the vacuum are given by the formula 

(3.9) U(w,z)= / x°( w ~ s > z ~ s ) "0( s ) ds, 



valid in each region p < 1 and p > 1 (as they avoid the point mass at w = z = 0), 
where ip : M — > R is an arbitrary, integrable function. A similar formula hold for 
the entropy flux. 

We recall that entropy solutions to the non-relativistic Euler equations admit 
bounded Riemann invariants, and therefore satisfy the restriction \v ± lnp| < C 
in the physical variables. In particular, the argument of the function / in the 
definition of y° remains in a compact set. Note also that the entropy kernel satisfies 

\X°(P,V)\<P 1/2 - 

3.3. Entropy kernel. We now return to the general model with e ^ 0. One of the 
main results in the present paper is the following characterization of the entropy 
kernel. 

Theorem 3.2 (Entropy kernel of the relativistic Euler equations). Consider the 
Euler equations for isothermal relativistic fluids (|2.1[) . The function 



fft(0,0;w,z), wz<0, 
10, «jz>0 ; 

is a fundamental solution of the entropy equation 

Xw Z + b(w + z) Xw + a{w + z) x z = -25 w=0 (g> S z=0 . 

It is solely a function of bounded variation and 



(3.10) 



lim x(w,z) = v " — -^—e 

Z—*0,WZ<0 Z 



l+e+{l-e)e e _ (1 _ e) =. w/4 



(3.11) 



ii m< v(^.) = 1 - £ + ( 2 1 + £)e£ % -^^, 

W— »U,lO2<0 Z 

Xw(w, z) = -{sgnz)- — £ + ^ + £ ' )e e- {1+e)2z/4 S w=0 + d{w, z) t wz<0 , 
Xz (w, z) = -(sgn w) 1 + £ + { l - £)e6W e^^ 4 5 Z=0 + C 2 (w, z) l wz<0 , 
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where C\^Ci are smooth functions. In the physical variables p,v, the kernel x = 
v) has compact support in the variable v G (— 1/e, 1/e) (for every fixed value 
p), and satisfies 

lim x(p, •) = 0, lim Xp(p> ') — ± _r a $v=o in the weak sense, 
p—>i p— >i± + e 

|X(P,U)I^P Q , a : = 2(1+2") ■ 

In the course of the proof of this theorem we will also show 

1 - e 2 1 + £ + (1 - e)e 
1 

and 



(3-12) ., 



4 / 2 

3 - 13 , , X 

1-e 2 / 1-e 2 \ l + e + (1-e e' 



C a («,,0) = 1 + ^w iZlZStZll e(1 -^W4 



Observe that the kernel x is only Holder continuous at the vacuum; the Holder 
exponent converges to 1/2 as e — > 0, allowing us to recover the behavior in ^J~p 
known in the non-relativistic regime. (In fact, as e — > 0, the expansion of the 
entropy kernel converges term by term to the one of the non-relativistic case.) 
Next, we introduce: 

Definition 3.3 (Notion of weak entropy). A pair of continuous maps (It, 5") : R+ x 
(— 1/e, 1/e) — ► R 2 is called a weak entropy pair of the isothermal relativistic Euler 
equations if the partial differential equations (13. ip hold in the sense of distributions 
and, moreover, 11 vanishes on the vacuum within any tame region, in the following 
sense: for any M > and e G (0, 1) there exist some constant C > such that 

(3.14) \U{p,v)\ + \2{p,v)\ <Cp a , (p,v)e7 e (M). 

For instance, the conservative and flux variables in (|1.1|) do qualify as weak 
entropies. From Theorem 13.21 we deduce: 



Corollary 3.4 (Entropy pairs of the relativistic Euler equations). In each of the 
regions p < 1 and p > 1 the formula 

(3.15) U(w, z) = / x( w ~ s i z ~ s ) V>( s ) ds 

Jw 

determines a family of weak entropies, where tj) : R — * R is an arbitrary integrable 
function. 

3.4. Entropy flux. Given an entropy U (with sufficient decay) we deduce from 
(I3.1[) , that the associated entropy flux is 



3(w, z) = / M{w', z) lL w (w / , z) dw' 

= X 2 (w, z) VL(w, z) — / \2 W (w' , z)U(w' , z) dw' . 
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A similar formula hold with A2 replaced by Ai and, therefore, 

&{ w i z ) =2 z ) + ^2{w, z)) U(w, z) 

1 f w 1 f°° 

— - / A2 W (ty ,«)lt(«; ,z)dw + - / A iz (w, z) U(w, z) dz . 

^ J — oo ^ J z 

These identities allow us to define the entropy flux kernel, denoted below by a — 
cr(p, v, s) from the entropy kernel x, and in turn to compute the entropy flux via 
the formula 

(3.16) 3"(p, v) = J(w, z) = / u{w — s,z — s)tp(s)ds. 



More precisely, the condition (J3TTJ) reads 

(3.17) a w = X2XW, <?z = Ai% z . 

By taking advantage of the Lorentz invariance property of the relativistic Eulcr 
equations we can decompose the entropy flux as follows: 

(3.18) a(p, v, s) =: u(v) x(p> v, s) + o-*(p, v, s), 

or, equivalently, cr(p, v, s) — ^1 x{p, v, s) + cr^{p, v, s). The key property (which is 
easily checked from the equations defining the entropy and entropy flux kernel) is 
that the dependence of x an d c with respect to s can be suppressed, since 

X(p,v,s) = x(w - s,z - s), 

cr"(/9, v, s) = <t"(iu — s,Z — s). 

Clearly, the condition (]3 . 1 7[) is equivalent to 

u w + z 1 
°iu = 7T~ Xw - XX + A 2 Xw, 

(3.19) _f / 

« w + z 1 
<*l = — Xz - ^X + Ai Xz 

Moreover we impose the following boundary condition on a* = a^(w, z): 

a* (0,0) = 0. 

In particular, the traces along the boundaries w = and 2 = can be determined 
explicitly from the boundary conditions for the Riemann function (|3.5p and the 
expansion l|3.11[) . By plugging (|2.4[) and (|3. . we obtain the following corollary 
of Theorem [ 



Lemma 3.5 (Entropy flux kernel). 

vl =C 5 (z) e -( 1+£ ) 2z / 4 S w=0 + C 3 (w, z)l wz<0 , 
al^C 6 (w)e^-^ w /U z= o + C 4 (w,z)t wz<Q , 

where 

C 3 (w, z) = (X 2 (w, z) - - - Z ) d(w, z) - -x{w, z), 
C 4 (w,z) = (Ai(w,z) - - ^ Z ) C 2 {w,z) - ~x(u>,z), 
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, l-e+(l + s)e ez . , l-e-(l + e)e e 
C 5 (z) := \z\ \ '— + (sgnz) 



4 v ° ' 2e 

_ . . . | 1 + e + (1 - e)eF w , , . l+e-(l-e)e E 
C 6 (w) := |u)| '- + (sgnw) ■ 



4 v o / 2e 

Finally, it is not difficult to check that, again in each of the regions p < 1 and 
p > 1, the formula 



^(uj,^) = / a(w — s, z — s)tl)(s) ds 

(3.21) 



(•u(«)x(it; — s, jr — s) + cr"(iy — s, z — s)) ^(s) oLs 
determines the entropy flux associated with (|3.15|) . 



4. Technical estimates 



In this section we give a proof of Theorem 13.21 The existence of the Ricmann 
function is standard, as it is determined by a Goursat problem for a linear hyperbolic 
equation with regular coefficients. The main issues to be dealt with are the behavior 
of x near p = and near v = ± 1/s, as well as the expansion of x along the boundary 
of its support. 

Lemma 4.1. The traces of the kernel along its support (wz < 0) are given by 



(4.1) x(w,z) 



B(z) := i-E+a+E)'" e-d+^V 4 , w ^0, 



and, moreover, within any compact set in p 

Proof. We need to integrate out the boundary conditions arising in (|3.5p . The 
second differential equation in (|3.5p implies 

1 — e 2 / /""'' 1 \ 

ln|3l(ti>',«;t«,«)H-C , o = — {-(w' + z) + 2 / ~-dv/) 

4 V J 1 — e z v / 

for some constant Co, where 

1 , f w ' 1 

dw' = / , . . — ; , . . dw 



1 — e 2 v J 1 — etanh (e(w' + z)/2) 

2 



In 



1-e 2 

The third equation in (|3.5p implies 

1-e 2 



((1 - £ )e £ ^'+^ + 1 + e) + + z )- 



ln|^( W ,0';w;,z)| +C = — —(-{w + z') + 2 / s-da/) 

4 V J 1 + e^i> / 

for some constant Co, where 

1 f z ' 1 



l + e 2 v J 1 + etanh (e(w + z')/2) 
2 



1-e 2 



In 



((1 + e)e< w+z '^ + 1 - e) + + z'). 
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Since Ji(w, z\ w, z) — 1, we find 
H{w l , z ; w, z) 

,-(l-e) 2 ( W '-w)/4 l- e +(l+ e ) e -(-» + '> =■ A(w' ■ W z' = Z 

(l-e)e=(™'-™)+(l+ e )e-^™+-) -A\W,W,Z), Z Z, 

e (l+e) 2 {z'-z)/i l+ e +(l- e ) e -«(" + ») f , v , 

^ (i +£ ) e e(,'-») + (i_ £ ) e - e (™ + ») W —W, 

and in particular 



#(0,0; 



e (l-e)V4 = A(0; W;0 ) = A(W), Z = 0, 

e -(l+efz/4 1- E +(1+ E )e" = ff( Q;0;Z ) = w = 0j 



which provides the desired behavior on \ along its support. 

Next, we write along the boundary z — (corresponding to u = R and w = 2R) 

A(w) = -p 5 ^ (1 + S + (1 - ^p 1 ^), 

which shows that x{Pj v ) i5 P° '■ O n the boundary w = we have it = — i? and 
z = — 2R and we find a similar estimate x(p> v ) < P Q since 

_B(z) = - j o 2 < 1 +- 2 ) (l - e + (1 + e)p 1+? 3 "). 

Since x is smooth, the estimate remains valid in any compact region. □ 

Lemma 4.2. The entropy kernel x{ w i z ) = #(0, 0; w, z)1 wz< q is a fundamental 
solution of the entropy equation 

Xwz + b(w + z)xw + a{w + z)xz = -2S W=0 ® <5 z =o- 

The traces of the derivatives Xw an d Xz along the boundaries z — and w = 
( while keeping wz < 0) are 

lim Xw{w, z) = Ci(w,0) = A w (w) 

z^O 

(42) l-£ 2 l- £ +(l+e)e E "' {1 _ Efw/4 _. ... . 

= — ] o e =-W(" , ) ! 



and 



lim Xz(w,z) = C 2 (0, z) = B z (z) 

w—>0 



^ _l^l + g+ (l-^ e , MV4 = _ l{z)B[zy 

Proof. In view of (|4. 1[) . we obtain 

X«,(w,z) = # w (0,0; w,z)l wz<0 - (sgn2;)B(2)(5 to= o 

= -(sgnz) l - e + ( l + e >" e -W''*8 m=0 + d(w, z)t wz<0 , 

(4.4) 1 

Xz(w, z) = !R Z (0, 0; w, z) l wz<0 - (sgnw) A(w) <5 2=0 

= -(sgnii;) 1 + £ + ( ^ £)e£ "' e^-^^/^^o + C 2 (w, z) t wz<0 . 
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On the other hand, from (14 . 1 j) we can determine 

Ci(w,0) = A w (w) 

= d , i + £ +(l-e)e™ (1 _ £) , w/i 
dw\ 2 

l- e +(l + e )e" e(1 _ sfw/4 = _- {w)A(w) 
8 

and 

C 2 (0,z)=B z (z) 

d fl -e + (l + e)e ez e _ (1+g) 2 z/4 



dz 



_ £2) l + e+(l- £ )e" e _ (1+e) ^ /4 = 



We differentiate (|4.4[) once again, and in view of (|4. 2[) and (|4.3|) we obtain 

Xwz(w, z) = R wz (0, 0; to, z) t wz<0 - (sgn to) R w (0, 0; to, 0) 5 z=a 

- (sgnz) B z (z) 5 W=0 -2 5(0)^=0 ® <5 Z=0 
= 3^(0, 0; to, z) l„, z<0 - (sgn to) A w (to) <5 Z=0 

- (sgnz) B z (z) S w=0 - 25 w = <S> <5 Z=0 

= ^ 2 (0, 0; to, z) t wz< v + (sgn to) a(to) A(to) 6 Z=0 
+ (sgn z) b(z) B(z) 5 w= q - 2 5 W = <g> <5 Z=0 , 
which allows us to compute 

X^z + + z ) Xw + a(w + z) Xz 
= $. wz (0,0;w, z)t wz<0 + (sgn to) a(w) A(w) 5 Z = 

+ (sgn z) b(z) B(z) 5 W = - 2 S w=0 <g> S z=0 

+ b(w + z) (31^(0, 0; to, z) \ wz< q - (sgnz) B(z)5 w = Q ) 

+ a(w + z) (3l z (0,0;w,z)t wz<0 - (sgnto)A(to) S z = ) 
= (3^u)z(0, 0; to, z) + b(w + z) 3l w (0, 0; to, z) + a(to + z) 3? z (0, 0; to, z)) t wz< o 

+ (sgn to) (a(to) — a(w + z)) A(w) <5 z= o 

+ (sgn z) (b(z) - b(w + z)) B(z) S w=0 - 2 6 W=Q ® S Z=Q 
= -2S w =o &> d~ z=0 . 

□ 

Lemma 4.3. T/ie traces of the derivatives Xw and \ z along the boundaries to = 
and z = ( while keeping wz < 0) are 

lim Xw(w,z) = Ci(0, z) 

1 £^ 2 , ,,x 1 - e + (1 + e)e ez _ n+F) 2 z/4 
(4.5) = — (1 - (1 - e 2 )^/4) 7j e z ' 4 

l-e 2 



(l-(l-e 2 )z/4) X (0,z) 
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and 

\imXz(w,z) = C 2 (w,0) 

z^O 

(4.6) = — (1 + (1 - e 2 )w/4) y - eS 1 £ > W ' A 

= (1 + (1 - £ 2 )^/ 4 ) X(V>, 0)- 

Proof. We give the proof for the boundary z — with u> > 0, the calculation for 
the other boundary being similar. Using (|3.6[) and the expressions (|3.4[) we write 

■j-Xz(w,0-) +a(w)xz(w,0-) = -b(w)xw(w,0-), 
aw 

where the coefficients a and b can be regarded as functions of w + z = w and the 
right-hand side is already known. By integrating the above equation we obtain 

Xz(w,0-) 

(4 7) r w - ' — — 

1 b(w') X w(w',0)e^ a{w " )dw "dw' + X z{V^)e- S ° a(w ' )dw \ 

where x 2 (0, 0) stands for lim„,^o+ Xz{w, 0— ). 

In view of gj}, x z (0,0) = -(1 - e 2 )/4. In view of g2J) and A(0) = 1, 

e-Jo"^']*' = A(w) = i±£±ii^£lf!!l e (i^) 2 »/4 

and 

J»'oK')*»" = A M = l + g+(l-e)e £W e _ (1 _ e) 2 K _ w)/4 
A(io') 1 + e + (1 - e)e ew ' 
Returning to (|4.7[) . we find 

6K)X-K,0)e/™'^")^'" = -6K)sK)^K)4rTr 



and 



thus 



= (1 g) (1 + e+ (1 - e)e™) e d-=) a ™/* 
32 

x»(w,o-) = (1 + (i - ^ 2 W4) AH, 

c 2 (,,o) = (i + (i - £ >/4) l±£±^£lf!: 



□ 



The singularities p. lip of the derivatives of x with its values at the boundary 
([312]) . ([3~13|) follow then from (gUJ) and from (|43|) . (|4~3f , (l4~5j) . and Moreover 
substituting (|3 . 11 [) into 

X P = (d w - d z )xR P , 

we obtain: 

Lemma 4.4. TVie traces of the derivatives Xp along p — 1± are 

!im i x P (p,v) = ±——2^=0. 
p-»i± 1 + e z 
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Lemma 4.5 (Derivatives of the entropy kernel and the entropy- flux kernel). We 

have 

(4 8) Xs ^ p ' v ' s ^ 

= X (1) (w - z) (-6 S = W + 6 S=Z ) + X ( - 2) {s -w,s- z) t( s - w )( s - z ) <0 , 

where 

X (> (w-z):=sgn(w-z) — - — e { e > { >' 

= sgn(w-z) lim x(€, v)k=w-z, 

and 

X i2) (p,v,s) = X {2 \s~ w,s -z):=- (d(w- s,z- s) + C 2 (w- s,z- s)). 
We have also 

(4.9) Js{P,V,S) 

= X (3) (w - z) {8 S=W + S s=z ) + X {i) (s -w,s- z)l( s _„;)( s - z )<o, 

where 

'\w-z\ l + s+(l~e)e^ w -^ 



X^(w-z) :=-0- 



2 2 



+S gn( W -,) 1 + £ -^-^' Z) ) e ^-»* 
= sgn(w-z)(\i(w-z,0)-—-—) lim y(f> rj)k=w-* 

and 

X {A) {p,v,s) = X {4) {s- w,s -z):=- (C 3 (w- s,z- s) + C 4 (w - s,z- s)) . 

Proof. In view of (|3.11|) and (|3.20p . we obtain 
d 

Xs(p,v,s) =— (x(w - s,z- s)) = ~(Xw + Xz)(w - s,z- s) 



Ci(w - s,z - s) l( w _ s )( z _ s )<o 

+ sgn(w-s) — e { e > y S)/ *5 Z 

C 2 (w - s,z- s) l(„,_ s )( z _ s )<0 



and thus 

Xs( P , v, s) =sgn (w - z) 1 + £+(1 2 ' £)e£KZ> e^ 2 ^ 4 (S s=w + S s=z ) 

- (Ci{w - s,z-s) + C 2 (w - s,z- s)) l( w _ s )( z _ s ) <0 - 

Defining XW( w — z) and X^ 2 \s — w , s — z) as in the lemma we obtain the derivative 
of the entropy kernel (|4.8[) . In view of (|3. 10|1 . we have also 

X (1) (w - z) = sgn(w - z) lim x(£> v)\i=w-z- 
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Similarly, for the entropy flux kernel we use using (13.201) and obtain 
d 

o*{p,v,s) = (a*(w -s,z-sj) = -(oi +(?t){w -s,z-s) 

l-£ + (l+£)e £ < z - s > , . l-£- (l + e)e £ ( z - s ) 
\z - s\ sgn (z - s) — 

-(l+e?(zs)/A _ Cz{w _ gjZ _ s )t {w _ s){z _ s)<0 

.l + e+(l-e)e< w - s ) , . 1 + e - (1 - e)e e( - w -^ 
\w-s\ sgn(w-s) — 

\w- z\ 1 + £+ (1 - e )e< w -^ . .l + e-{l-e)e< w - z "> 
— 2 + Sgn {W ~ Z) 2e 

e (i-^-^ {Ss=w+Ss=z ) 

- {C 3 {W - S, Z - S) + C 4 (w -S,Z- s)) l( w _ s )( z _ s )< - 

Defining X^ (w — z) and X^ (s — w, s — z) as in the lemma we obtain the derivative 
of the entropy- flux kernel (|4.9|) . Furthermore, in view of (|2.4p . (|3.10p 



X^ =sgn(w-z)(\ 1 (w-z,0)-— -) Urn x(£,v)k=v>-» 

Z J)^0,(;r/<0 

= sgn(-u) + z)(A 2 (-w + z,0) ) e (1 e) ( ™ z)/4 . 

□ 

5. Uniform estimates for the Lax-Friedrichs scheme 

In this section, we follow DiPerna (1983) who considered non-relativistic poly- 
tropic fluids, and we apply the Lax-Friedrichs scheme. As observed by Hsu, Lin, 
and Makino (2004), DiPerna's arguments carry over to relativistic fluids. Our main 
purpose is to derive new uniform bounds, and establish that Lax-Friedrichs approxi- 
mations remain in a tame region, which is uniquely determined from the initial data 
of the initial value problem under consideration. 

We begin with: 

Lemma 5.1 (Riemann problem). Consider the relativistic Euler equations (jl.lj) 
for e G (0, 1]. The Riemann problem corresponding to an initial data made of a sin- 
gle jump discontinuity (at x = 0) separating two constant states (pi,vi) and (p r ,v r ) 
admits a unique self-similar solution (p,v) = (p,v)(x/t) satisfying all entropy in- 
equalities. The solution of the Riemann problem satisfies the uniform L°° bounds 
in the Riemann invariant variables, 

< W(x/t) < W := sup(Wj, W r ), 

(5 ' X) < Z{x/t) < Z Q := sup(Zj, Z r ), 

where the notation W, Z stands for the modified Riemann invariants associated with 
the variables (p,v). As a consequence, there exists a constant M > depending 
only on (pi,vi) and (p r ,v r ) such that 

< p{x/t) < M, 
l-e\v(x/t)\>(p(x/t)/My', 
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where e' was defined in the introduction. 

Proof. When e S (0, 1) the Riemann problem was solved by Smoller and Temple 
(1993). In particular, they established that the shock curves lie in the interior of the 
rectangular regions limited by the integral curves of the characteristic fields. This 
property implies directly the maximum principle (|5.1[) . Recall that the rarefaction 
curves are given by 

Ri : u + R = in + i?;, P < Ph 
R 2 :u- R = ui- Ri, P>Pi- 

Hence on Ri, w is constant and z increases and on R 2 , z is constant and w increases. 
The shock curves are given by (i = 1, 2) 

Si : p = p(v; Pl , Vl ) =pi+pi (3{v, «,) (l + (-l)Vl + 2//?) , 

where 

_ (1+e 2 ) 2 (v-vtf 
Pl ' ij ' 2 (l-eV»)(l -£*«?)• 
From the definition (12.51). the condition (15.11) can be rewritten in the form 



^ P m^ < (i±£^) 1/(2e) < wop( x/ t)-y^\ 

Zo V 1 — sv(x/t)/ 

l + £ 2 

which leads to (|5.2|) with a constant M s := (sup(Wo, Zq)) depending upon e. 
In fact, from p = {W Zf 1 ^/ 2 , 

p(x/t) < (W ZoY 1+£2)/2 < M E . 

On one hand, [t^J/l] < p{x/t)- 2 ^( 1+ ^W 2e leads to 

2 2 

1 ~ £V{x/t) ~ l + p(x/t)-^(^)W ^ ~ l + (p( x /t)lM e )-^^ 
> (p(x/t)/M e ) £ ' 
On the other hand, the inequality 

l + ev[x/t) ^ „^ M 2 £ /(l+e 2 ) 7 -2e 



1 — ev(x/t) 

leads us to 



> p(x/ty e '^ +e >z 



1 + ev{x/t) > — — — —■ —z - > 



1 + p(x/t)-W+*)Z? ~ i + { p ( x / t )/M e )- 2£l{1+ ^ 
> (p(x/t)/M e f . 



This establishes (|5T2|) . 

Now, by expressing Wi, Zi,W r , Z r in terms of the initial density and velocity, 
one see that the constant M £ can be taken to be independent of e. In fact, from 
the definition (12.51) . 

l/(2<0 . / -i r \ -l/(2e) 

V 1 - e«o 



V 1 - ev J 
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SO 



/ii \ (l+e 2 )/(2e) 

\l-ev J 



, A 1 + £ , / 1 + ev 



Poe with A = — - — log 



2s \ 1 — evo 
By taking ^(z) = £ log (±±§ ) (x ^ 0) with p(0) = 0, 



^ = (1 + e 2 )v a ip(ev Q ) 
= (l+s 2 )v { 
and p'(0) = 0, ip"(0) = 8/3. Hence 



(l+e 2 )v a ^^(ev ) 2 (l + o(\ev \)) 



= Po cxp ^(1 + £ 2 )« ^(e«o) 2 (1 + o(|^oD) 

Similarly, we can show Zq +£2 = /9 exp ((1 + e 2 )v ^^-(ev ) 2 (1 + o(|et>o|))) ■ There- 
fore M e can be expressed only in terms of (pi,vi) and (p r ,v r ). □ 

Consider a family of cartesian discretizations of the spacetime K+ x R, based on a 
time length t and a space length h, where the ration r/h is kept fixed while h — > 0. 
Set i n := nr (n being a positive integer) and xj := j h (j being an integer). The Lax- 
Friedrichs scheme allows us to construct approximate solutions 
to the relativistic Euler equations, which, for any pair (n,j) with n + j even, is 
constant equal to p™,Vj in every slab [t n ,t n +i) x Ij := [t n ,t n+ i) x (xj-i,Xj + i). 
The initial data po,v is averaged over each initial cell 

{p%v°A := -!- / (p Q ,v Q )(x)dx, j even. 



Given a piecewise constant approximation at a given time t n , we solve a Riemann 
problem in the neighborhood of each point Xj + i, for j such that n + j is even, and 
we then average the solution at time t = t n+ \ — over the intervals (xj 7 Xj +2 ): 

(P]tl v jti) := ^hl ^ ^XWi - °> x ) dx ' n + 3 even - 

Observe that, if e e (0,1) (as well as for e = 1), the characteristic speeds 
Ai, A2 remain bounded globally, even when the velocity approaches the light speed. 
Indeed, all wave speeds in the problem under consideration are bounded. (Note that 
a quite different situation is met with the non-relativistic model corresponding to 
£ = 0, for which the characteristic speeds are unbounded.) To avoid any interaction 
in two neighboring Riemann problems, it is necessary to restrict the ratio r/h by 
the CFL stability condition 

/ x T ( \» h \ +1 \ 1 
(5-3) — max ^— =- < -, 

where the maximum is taken over all (t, x). For instance, a sufficient condition is 

<"» 51 < 5- 

which clearly becomes more restrictive as £ approaches 0. 
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We observe that both the Riemann problem and the projection step satisfy 
uniform stability estimates. Finally, from Lemma 15. II we deduce the key estimate 
of this section: 

Lemma 5.2 (Uniform a priori bounds). Provided that the CFL condition (|5.3|) 
holds, the approximations (p h ,v h ) satisfy the tame condition 

< p h (t,x) < M, 

(5.5) 

l-e\v h {t,x)\ > (p h (t,x)/M) £ , 

for some constant M > depending solely on the initial data. 

In addition to the uniform amplitude bound on p , v h one can also derive entropy 
dissipation bounds which imply that the entropy inequalities associated with any 
weak entropy are satisfied up to an error term vanishing in the distribution sense as 
h — > 0. We omit the details and refer to DiPcrna (1983) and Hsu, Lin, and Makino 
(2004) for further details on Lax-Friedrichs approximations. 

6. Reduction of the Young measure 

We rely on the theory of compensated compactness for nonlinear hyperbolic 
systems for which we refer to Tartar (1979 k 1983), Murat (1978 & 1981), and 
DiPcrna (1983). In view of the uniform bounds derived in Lemma 15.21 we can 
associate to the sequence p h ,v h a Young measure v = v t . x supported in a tame 
region T(M) for some uniform constant M > 0. By definition, for almost every 
point, v t .x is a probability measure in the variable p,v which has compact support 
in the Riemann invariant variables W, Z. The Young measure allows us to compute 
the weak limit of any composite function of the sequence p , v , that is: 

(6-1) g(p\v h )^(v,g) 

in the sense of distributions, for every function g that is continuous in the Riemann 
invariants W, Z. 

Relying on standard arguments one can check that the entropy dissipation mea- 
sure associated with a weak entropy pair (U, 9"), 

d t U(p\v h )+d x 3(p h ,v h ) 

belongs to a compact set of the Sobolev space H^. By the div-curl lemma, this 
property implies that v satisfies Tartar's commutation relation 

(6.2) (v, U^ 2 - U 2 5i) = (v, Ui){v, J 2 ) - (v, U 2 ) (i/, Ji) 

for any two weak entropy pairs (Hi, 5i), (U2, 3^2 )• Plugging the entropy-entropy flux 
pairs given by (|3.15|) and (|3.21|) in the relation (|6.2j) and dropping the test-function 
■0, we obtain for almost every (t, x) and for all s, s' 6 R 

(6.3) W,x(s)a(s') -x(s)o-(s')) = (v, X (s))(v,o-(s')) - (u, X (*')) 

Equivalently, this identity holds with a replaced by cr" = a — u \- 

Our main result in this section is as follows. We fix a point (t, x) (where the 
above relation holds) and consider the Young measure at that point. 
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Theorem 6.1 (Young measure reduction for the relativistic Euler equations). A 
probability measure v with compact support in the (W, Z) -plane and satisfying (|6.3j) 
for all s, s' is either a Dirac mass 



In other words, there exists a function (p, v) = (p,v)(t,x) defined for almost 
every (t, x) such that, almost everywhere, either p(t,x) > and v t x = o~( p , v )(t,x) or 
else p(t, x) = and v ttX is supported in the vacuum line (the velocity begin then 
irrelevant). This immediately implies, for instance, that p h as well as p h v h converge 
to their limits p and pv, respectively. By a (standard) property of consistency of 
the Lax-Friedrichs scheme the limit must be a weak solution of the relativistic Euler 
equations satisfying entropy inequalities for all weak entropy pairs. This completes 
the proof of our main result stated earlier in Theorem 11.21 

Proof. We closely follow the method used by LeFloch and Shelukhin (2005). Given 
some si € M, it will be convenient to write \i = x(p> v i s i)- Given si,S2,S3 S K, 
we consider the corresponding identity (|6.3j) for the three pairs 



We then multiply these identities by {v, {v, Xi)> { v i X2), respectively, where Xj — 
x(p, v, Sj) (j = 1,2,3), and we add up these three identities. Then, due to the 
symmetry of the expressions, the sum of the right-hand sides vanishes identically: 



{v, X3) {{v, Xi) {v, o\) - {v, X2) {v, o\)) + {v, xi) {{v, X2) {v, o\) ~ {v, X3) {v, cr|» 
+ W, X2) ((v, X3) crj) - (f, Xl) (v, (T3)) = 0, 



(f,X3)(^Xl CT 2 - X20-\) + {v,Xl){X2°\ - X3°~ 2 ) + (v,X2)(v,X3<Tl - Xl^l) = °- 

Then, by differentiating once in S2 and in S3 and by setting dj := d/dsj, we 
obtain 

{v,dzXs){v,Xifo<A - foxzoi) + {v,d 2 X2)(v,d 3 X30-l - Xi9 3 <rf ) 

(6.4) 

= - (f,Xl)<32X2§3<73 ~ d 3X3d 2 Crl), 

which is valid in the sense of distributions in M 3 . In view of (|4.8[) and (|4.9p . by 
setting :— X^\ s=Sj we find, on one hand, 



v = 5, 



or else has its support included in the vacuum line 



supp v C {p = 0} = {WZ = 0}. 



{si,s 2 ),(s 2 ,s 3 ),(s 3 ,s 1 ). 



where aj — cr^(p, v,Sj), whereas the sum of the left-hand side is 



(6.5) 



Xid 2 <T 2 - d 2 X20-{ 




where Ej := {(w — Sj)(z 




(6.6) 



<9 3 X30i - Xid 3 03 
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and, on the other hand, 
(6.7) 

#3X39202 - #2X2 d 3 a\ 
= 2X^X<- 3 \S S2=W 5 S3=Z - 6 S2=Z 5 S3=W ) 

+ ((x^x^ + xWxjP) s S2=w + (x^x^ xVxjP) S S2=Z ) t E3 

+ (( - X^X™ - X^X™) 5 S3=W + (X^X^ X^X^) 5 S3=Z ) 1 

+ (x^x^-xPx^)t E2 t E3 . 



In view of the formulas (|6.5[) and (|6.6p . the right-hand side of (|6.4|) contains 
products of functions with bounded variation and Dirac masses, plus regular terms. 
Such products were already discussed by Dal Maso, LeFloch, and Murat (1995). 
On the other hand, the right-hand side of (|6.4p is more singular and involves also 
products of measures. Our calculations will show that the left-hand side of (|6.4p 
tends to zero in the sense of distributions when S2, S3 — > Si, while the right-hand 
side tends to a non-trivial limit. 

We test the identity (|6.4|) with the following function of the variables si, s 2 , S3 



(6.8) i/tfa) <p%(si - s 2 )ip1{s 1 - s 3 ) := ip(sx) -^^((si - s 2 )/k) <£ 3 ((si - s 3 )//e), 

where k is a small parameter, tp is a smooth and compactly supported function, 
and ifij : R — * R (j = 2, 3) are mollifiers such that 



^ Vj(sj) rfsj = 1, supp <pj C (-1, 1). 



W > 0, 

We consider first the right-hand side of (|6.4p . Applying the test- function ([6 
to the distribution {y,Xi){^2X3d 3 cr 3 — d 3 x 3 d 2 <r 2 ), we obtain the integral term 

(6.9) / (^, Yi)(<9 2 X23 3 03 - d 3 X3d 2 al)ip(s 1 )tp 2 i (s 1 - s 2 )(f 3 (s 1 ~ s 3 ) dsids 2 ds 3 , 
it 3 

4 

which we decompose as a sum ^ Relying here on Lemma 14.51 we distinguish 

i=l 

between products of Dirac measures 

Jf := 2 f ^(«i)<i/, xi) (l (1) I (3) (^(si - t«)^(«i - 2) - ^( Sl - zM( Sl - «;))) d fll 
products of Dirac measures by functions with bounded variation 

: JR K ' 

</?3( s i _ s 3 )\ Ez ds 3 ) dsi 

tP( Si )(v,Xi)(<P k 3 (si-w) I (X^X^+X^X^) 

l JR K ' 

<p 2 {si - s 2 )t E2 ds 2 ) dsi 

— ■ 1 2A — X 2,2: 



1 2 
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/;; = / i,( Sl ){v, X i)(rt(si-z) f (xWxW-xWxW) 

^(si - s 3 )t Ea ds 3 ) ds! 

^MfaxiXtfO*-*) / (x^x^-x^x^) 

^2 ( s i - s 2 )t E2 ds 2 ) dsi. 

The remainder 

11 := f ^( 8l )<i/,xi>< (X^X^ X< 2) *f) t E2 t E3 ) 

<P 2 {si - s 2)^3(si - s 3 )dsids 2 ds 3 

involve functions only. 

By a change of variable we see that the integral term 

Ii=- [ [ [ iP( K y + w) X (p,v, K y + w)X^X ( V 

K JW,Z JW',Z' Js. 

{vZivMiv + {w- z)/k) - tp$(y + {w- z)/k)<p$(jj)) dydvdv 1 
tends to zero, i.e. If — > 0. The same is true for the smoothest term I£, in fact 



( X^(p,v lS3 )^ 3 ( Sl - S3 )ds 3 xW(p, v,s 2 ) 



W,Z JW',Z' 

<^(si -s 2 )ds 2 - I X {2 \p 1 v,s 2 )Lp 2 l {s l -s 2 )ds 2 I X^(p,v,s 3 ) 



T K 



Vzisi - s 3 ) ds 3 jip(si)x(p,v, si) dsidudv' , 
which clearly tends to 

f f f (xWfoVssjxWip^Si) - xW(p,v, Sl )xW(p,v, Sl )) 

JW,Z JW,Z> JR K ' 

ip(si)x(p, v, s\) dsidvdv' = 0. 
Next, let us consider the term 1% = 1% 1 — Ig 2 in (|6.9p . By defining 

Q" :=X( 2 )X( 3 )+X«X( 4 ), 

<?7 := Q~(p> v > s /)' x' ■= x(p', v ', s ),x'j ■= X(p',v'sj), 

we can write 

I 21 = -[ [ [ ^(si)xi^(si - w) [ Q 3 f 3 (si- s 3 ) ds 3 d Sl dvdv' '. 

JW.ZJW'.Z'JR Jw 

Lemma 6.2. Let f,F:~EL—>-M. be continuous functions. Then, for every interval 
[a, b], [a', b'] C K, the integral term 

rb' rb 

I K (a,b,a',b') := / /(si)^(si - a) / F(s 3 )¥>£(«i - a 8 ) ds 3 dsi 

■/a' */ a 

converges (when k — > 0^) toward 



/(a) F(o) (A 2i3 1 a— a a = 



:6' 
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where the coefficients A 2 $ := B 2 . 3 + 6*2,3 and B 2i3 and 6*2,3 depend on the mollifying 
functions as follows: 



roc ryi 

B 2 , 3 := / <P2(yi)V3(y3)dy 3 dyi, 

JO J -00 

/0 ryi 
/ (p2(yi)(fi3(y3)dy 3 dyi. 
-OO J —OO 



Formally the integral I K has the limit 

I(a,b,a',b r ) :=/ /(si)5 8l=a / F(s 3 )J S3=Sl ds 3 dsi. 

J a' J a 

Proof. Making first the change of variables S3 = si — ny 3 and then s\ = ny\ + a, 
we can write 

n(b'-a)/n ryi — {b-a)/n 

I K = - f(nyi +a)ip 2 (yi) / F {n^ - y 3 ) + a) ip 3 (y 3 ) dy 3 dy 1 . 

J(a'—a)/K Jyi 

Clearly, we have / K — > when a < a' or a > b' . 
Now, if a = a', we can write 

p(b'-a)/K ryi— (b-a)/K 

I K = - / fi^yi + a)(p 2 (yi) F (n(yi - y 3 ) + a)) cp 3 (y 3 ) dy 3 dyi 

Jo J Vl 

poo pyi 

f(a)F(a) / <f2(vi) ^3(2/3)^3^2/1 

as k -* 0. The other values of a can be studied by the same argument and this 
completes the proof of Lemma 16.21 □ 

Applying Lemma HT2l we find that 1 tends to 

4>{w)x' {w)Q~ {w)(A 2 ,3^-oc<w<oc + B 2 . 3 1 W= „ 00 + C 2 . 3 t w=OQ ) dvdv' 

W,Z JW',Z' 

= -A 2 , 3 [ f i){w) X '(w)Q-(w)d V dv' 

Jw,Z Jw.Z' 

and that I 22 tends to 

-A 3t2 [ [ ^{w)x'{w)Q~{w)d V dv'. 
Jw,z Jw,z' 

We conclude that the limit of I 2 is 

A [ [ i){w)x'{w)Q-{w)dvdv' , 
Jw,z JW,Z' 

where 

A :=A 3i2 - A 2 . 3 

l-Vl 

(y>2 (2/3)^3(2/1) - <P2 (2/1 V3 (2/3)) dy 3 dy x . 
We can determine similarly that the limit of 1% which is found to be 
-A I f 4>(z)x'(z)Q + {z)dvdv', 

JW.Z JW'.Z' 
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where 

In conclusion, we have identified the limit of the term (16.91). as 



A[ f [ i/>(w)x' (w)Q~ (w) dvdv' - [ [ ^{z) X '(z)Q + (z) dvdv' 

\Jw,Z JW',Z' JW,Z JW'.Z' 

One can also check (LeFloch and Shelukhin, 2005) that the distribution 

{v, d 3 xa) (v, Xi<h<A - d 2 X2o\) + (u, 82X2) (f, 83X3^1 ~ Xi9 3 ^) 

(which is the left-hand side of (|6.4[) applied by the test-function (|6.8p tends to zero 
as k — ► 0. 

Since the molifying functions (f2 and if3 can always be chosen such that A ^ 0, 
we have reached the following conclusion 

tj}{w)x' {w)Q~ {w) dudv' - / / iP{z)x'{z)Q + {z)dvdv' 
(6.10) Jw,z Jw,z' Jw,z Jw'.z' 

= {v®v l ,ii{w)x'{w)Q~{w)- ^{z)x'(z)Q+(z)) =0. 

By Lemma 14.51 we have 

Q~(w) =Y (2) (0, w - z)X (3) (w -z) + X {1 \w - z)Y (4) (0, w ~ z), 

Q+(z) =X^ (-w + z, 0)X<- 3 \w -z)- - z)X& {-w + z, 0). 

Observe that the test-function tp is arbitrary and can be used to localize the 
equation 16 . 1 0|) . In turn, we end up with a necessary condition satisfied by the 
Young measure v: 

(v ® v', Eg) = 0. 

Here, v®v' denotes the tensor product of v and v' (another copy of the Young mea- 
sure u), while g is some (non- negative) characteristic function, while the function 
E is defined by 

E(p) := lim y(u?, z) (at(w, z) + a\(w, z)) + (xw(w, z) + Xz(w, z)) J(w, z), 

p and w being related by w = 2R = 2(mp)/(l + e 2 ). An analogous identity holds 
with lim replaced by lim . For convenience, we write simply x(w,0) 

instead of lim x( w i z )? e ^ c - This allows to rewrite the expression of S as 

z—>0.wz<0 

E{w) = - (xK0,0)<7*K0,0) + X s(w,0, 0)o*(vj,Q,0)) 

,() ' " ' = X {w, 0) (ai(w, 0) + 4(w, 0)) + (*„(«;, 0) + x ,{w, 0)) o*(w, 0). 

Provided the coefficient E keeps a constant sign, the above condition implies 
(LeFloch and Shelukhin, 2005) that the support of v either is a single point, or 
else has its support concentrated where W Z = 0, which is nothing but the vacuum 
line. Recall that the relativistic equations are automatically satisfied if the density 
vanishes identically. 

It remains to establish that E(p) < 0. More precisely, we only need this to 
hold for sufficiently small p, since by using the scaling invariance property of the 
relativistic Euler equations (Lemma ll.3|) we can always ensure that the range of 
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the Lax-Friedrichs approximations and therefore the support of the Young measure 
is included in a neighborhood of the vacuum. We will actually prove that 

(6.12) S(p) < -i/? 2Q < 0, 0<p<l. 

We set 3(w) =: lim x 2 {w , z) Vl(w) , and we observe that lim xi w i z ) 
behaves like p a near the vacuum and, therefore, we need to check that 

(6.13) «X-)-(-i + !)+<**> 

< 0. 

The term 0{e 2 ) should be bounded by a constant times e 2 , uniformly for all large 
(negative) w. This will show that f2 has a fixed sign for all sufficiently negative 
values of w, and this will indeed complete the proof of the theorem. 

It remains to determine the sign of the function 5. For the sake of comparison, 
we recall first the relevant formula when e = 0: 

A?K0) = |-1, A§(«;,0) = ~ + 1, 



and 



x°(w,o) = e w ^, 

X °Jw,0) = X z(w,0) = - J(l + w/A)e w '\ 

(x!+X°)K0) = -^e^ 



and 



<T>,0)=(|-l) e-'V't^O) 



_ e w/4 



2 16 32 
Hence, in the non-relativistic case we do have 

For the general case e^Owe have already determined the traces of the entropy 
kernel along the boundary z = while keeping wz < (see (|4.1[) . (|4.2p and Q4.6[l ): 

0) = (1 - e + (1 + £ )e^) e^ 2 ™/ 4 

= 1 ~ £ " 1 - £ + i 1 + £ \ eSW X {w, 0) ^ a(w) X (w, 0), 
4 l + e+(l-e)e £U ' v ; 1 Ml ' ; ' 

0) = (1 + (1 - e 2 )w/4) X (w, 0). 



(6.14) 
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On the other hand, for the entropy kernel, from (|3.17() . (|2.4|) and (|4.2j) it follows 
that 



a w (w,0) = X 2 (w,0)xw(w,0), 

= (l-e) 2 «./4 



1 1 - e - (1 + e)e<™ 1 - e 2 1 - e + (1 + 6)6™ 



e 1 - £ + (1 + e)e ew 4 2 
8s 



that is 

cr(u>,0) = / <r w (w' , 0)dw' 



_ g2) l- £ -(l+ £ ) e - e(1 _ £)v/4 du/ 

, 8£ 



So, we find 



g ( m ,0) = - 1 + E " (1 ~ E)e "e"^ 
v ' ; 2e 

11 + £ — (1 — £~)e £W 

= — il a. 1 — b^ w <°) = AiK0)xK0) 
e 1 + e + (1 — eje 6 ™ 

and thus in view of (|3.19|) 

(6-i5) ; x 

ai(w,0) = (X 2 (w,0) - |) x»K0) ~ 2 X(U; ' 0) - 
Furthermore, from ()3.19j) we get 

(6.16) a»( w ,0) = (Xi(w,0) - |) x ,(tu,0) - \x{w,G). 

Finally, we are in a position to compute the quantity of interest. In view of 
(|Q3|) . (11113, and (EUSl), we can obtain 

S W W ■■=x(w,0)<jI(w,0)+ Xw (w,0)<jHw,0), 
H«(p) := X ( W , 0) 4(w, 0) + 0) a\w, 0), 

as follows: 

sW(p) = x(u,,o) (Va^o) - 1) X ,(tu,0) - ^xKo)' 
+ x»Ko) (a!Ko)-|)xKo) 

= x («,,0) ^2 (Ai( W ,0) - I) X z(w,0) - ^x(<W 
= -\x 2 {w, 0) ((1 - e 2 ) (1 + (1 - e 2 W4) (a x («;, 0) - |) + l) 
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and 



3 W (P) = XKO) (^X 2 (w,0) - |) Xw (w,0) - i X K0) 
+ *«(«>, o) (A!K0)-|)x(w,0) 
= xK 0) f(Ai («;, 0) + X 2 (w, 0) - «;) 0) - \x{™, 0) 

= X (™,0) 2 ((\i(w,0) + X 2 (w,0) - w)a(w) - X - 



Thus, we conclude that H(io) = x 2 (w,0) fi(iu), where 

Sl(w) := - ^(1 - e 2 ) (1 + (1 - £ 2 )W4) (AiK 0) - |) 
+ (Ai(u>, 0) + A 2 (to, 0) - w)a{w) - 1. 



(6.17) 
We easily see that 



Ai(w,0) = w/2 - l + 0(e 2 ), 
\ 2 (w,0) = w/2 + l + 0(e 2 ), 
a{w) = -l/4 + 0(e 2 ), 
thus (|6.13p holds and S vanishes only at the vacuum p = 0. □ 

The formulas derived in the present paper converges formally to the ones of the 
non-relativistic case, as we now show by using the notation introduced by LeFloch 
and Shelukhin (2005). Just before the identity (5.7) on p. 424 of that paper, the 
function D{R) := Q~ (w) is defined while Q~{w) :— G x — G h is introduced on 
p. 420. The terms G x and G h are introduced in Theorem 4.6 on p. 414, that is: 

G x (R,u) = -2\R\f'(0)e R ^ 2 , u<-\R\, 
G h (R,u)=e R/2 (2R+^), \u\ > \R\. 

The calculation in the proof of Theorem 4.6 should be modified, as follows. 
First of all, given a test function (p = <p(s), we can write 

(x\ <P) = - X(s)<p'(s) ds = -e R / 2 / v'(s)f(\u - s\ 2 - R 2 ) ds 

JR Ju-\R\ 

= -e R /\y{u+\R\)- V >{u-\R\)) 

ru+\R 



2<p(s)f'(\u- s\ 2 - R 2 )(u-s)ds, 

u-\R\ 

which yields 

Xs = e R/2 (5 s = u -\r\ - 5 s =tt+|R|) + G x (R,u- s)l| u _ s |<|fl|, 
where, for all \v\ < \R\, 

G X {R, v) = -2e R/2 vf'(v 2 - R 2 ). 

Hence, we find 

G x (R,u) = 2\R\f(0)e R / 2 , u < -\R\, 
and this formula contains a plus sign, instead of a minus sign as stated originally. 
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Second, for the expression of the derivative of the entropy flux h s , a term f'(0) 
should be added, as follows. We write 

h s = sgn(w- s)e~ |,l ~ s|/2 ((5 s=tl _| i? ,| - S S=U+ \ R \ + ^sgn(u - s)l\ u _ sMR ^ 

~2 °( {u - s)e r/2 f{lu _ sl 2_ r 2 ) \ dr 

J -(\R\W\u-s\) ° s V 7 

+ 2/'(0)e-l«- s l /2 ( W - s )sgn (u - s) - 2 e ^-*l /2 l| tl _ s |>| i? |/'(0)| W - s\ 
= e R l 2 {5 S=U ^ R \ + 6 S = U+ \ R \) + -e-l"- s l /2 l|„_ s | < | i? | 



2 

' (e r ' 2 f'{\u - s\ 2 - r 2 ) + 2e r ' 2 \u - s\ 2 f"{\u - s\ 2 - r 2 )) dr 

-(|R|V|«-s|) V ' 

+ 2/'(0)e-l u - s l/ 2 | U - S |l| u _ s |<| fl | 



+ 2 



and, therefore, 

h s = e R ' 2 (S s=u _ m +5 s=u+lRl ) +G h (R,u- a)l| u _ s |<| fl | 

with 

G h (R,v) : = e-l"l/ 2 (i+2/'(0)M) 



Hence, we find 



e -^(i-iH). 



G>-tR,u) = e ""(l + f), u<-\R\ 



In conclusion, for small p we have R < and 

= e-/ 2 (-i + f)^e-/ 4 (-l + ^), 

which is precisely the same expression as the limit e — > of the expression (|6.13|) 
obtained in the present paper. 
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